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Abstract. We have analyz ed the catalog of 2106 non-kinematically selected metal poor stars in the solar neighborhood pub- 
lished bv lBeers etafl feoOOX with the goal of quantifying the amount of substructure in the motions of the fastest halo stars. We 
have computed the two-point velocity correlation function for a subsample of halo stars within 1-2 kpc of the Sun, and found 
statistical evidence of substructure, with a similar amplitude to that predicted by high resolution CDM simulations. The signal 
is due to a s mall kinematic group whose dynamical properties are compared to the stellar "stream", previously discovered by 
iHelmi et alj ll999l) . If real, this high velocity moving group would provide further support to the idea that substructures remain 
as fossils from the formation of the Galaxy as expected in the CDM scenario. 
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1. Introduction 

Observations of metal-poor stars have many times been 
used to discriminate (although not yet conclusively) 
bet ween alternative galaxy format i on scenarios (e. g ., 
lEggen. Lvnden-Bel l & Sandagel 1 19621 ISearle & Zirrnl ll978t> . 
In recent years, considerable effort has been put into under- 
standing the properties of galaxies within the hierarchical 
paradigm of structure formation in the Universe which, so far, 
looks to be the most successful theory. In a CDM Universe, 
the first objects to form are small galaxies which then merge 
and give rise to the larger scale structures we observe today. 
Thus, structure formation occurs in a 'bottom up' fashion. 
This theory predicts the presence of substructures (tidal tails, 
streams) due to the mergers and accretion that galaxies have 
experienced over their lifetime. 

Direct comparisons to observations have shown that this 
model can reproduce the properties of both the local and the 
distant Universe. Several examples of mergers and galaxy in- 
teractions have been observed in the Milky Way, such as the 
disrupted Sagittarius and Can is Major dwarf galaxies (e. g., 
Ilbata et alJll995t iMartin et alJl2004l) . the phase-space st ream 
of halo stars in the solar neighborhood ( Helmi et al. 1999), and 
the ring in the outer Galaxy INewberg et al J l 2002). Similar 
substructures have also been found in the halos of other nearby 
galaxies, such as M 31 JFerguson etai1l2002l) and NGC 5907 
( Zhe ng et al.l l 1999). showing that accretion may be a common 
phenomenon in the evolution of galaxies. 
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Although the stellar halo accounts for only about 1 % of the 
luminous mass, it plays a crucial role in studies of formation 
and evolution of the Galaxy. Signatures of the hierarchical na- 
ture of galaxy assembly are expected to be most obvious in 
this component. Moreover, stars in the halo are generally very 
old and metal poor, i.e. they can be considered more pristine. 
These are in fact the stars thought to have been formed in satel- 
lite galaxies that merged to form our Galactic halo (Robertson 
et al. 2005). 

At present, the best measurements of the halo kinematics 
are obtained from the analysis of samples of stars located in the 
solar neighborhood. Studies of the kinematics of various stellar 
populations in the Galaxy have long been limited - especially 
for the inner Halo - by the lack of large samples of stars with 
accurate distances, metallicity and kinematics. 

Beer s et alJ lEo OO) compiled an extensive catalog of metal- 
poor stars selected without kinematic bias, and with available 
proper motions, radial velocities, and distance estimates for 
stars with a wide range of metal abundances. In this paper we 
analyze this data-set, which has already provided support for 
constraining plausible scena rios for the formation and e volu- 
tion of the Milky Way (e. g., IChiba & Beerj2000lEool . 

The layout of this paper is as follows: In Sect. 2, we assem- 
ble a sample of metal-deficient ([Fe/H] < -1.5 dex) ha l o stars 
up to 2 kpc of the Sun selected from the iBeers et alJ (|2000) 
catalog. In Sect. 3, we explore their phase-space distribution 
and quantify clustering by means of the two-point correlation 
analysis. In Sect. 4, we compare to theoretical predictions from 
CDM simulations. In Sect. 5, we further examine evidence for 
kinematic substructure in the space of adiabatic invariants. The 
conclusions of the present study are given in Sect. 6. 
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2. Halo samples from Beers et al. (2000) catalog 

Beers et al. (2000) presented a revised large catalog of 2106 
metal poor stars in the solar neighborhood selected without 
kinematic bias. Within this sample, 1258 stars have distance es- 
timates and radial velocities as well as proper motions, so that 
the three components of the space velocities can be derived. 

Throughout this work we adopt Galactic velocity compo- 
nents relative to the Galactic center, where as usual U, V and 
W are positive toward the center, in the direction of the rota- 
tion, and toward the north Galactic pole, respectively. We as- 
sume that the local standard of rest (LSR) rotates with a ve- 
locity of 220 km s _1 about the Galactic center, and that the 
peculiar velocity of the S un relative to the LSR i s given by 
(10.00. 5.25. 7.17) km s" 1 JPehnen & Binne\ll998l) . 

2.1. Selection criteria and selected sets 

Since we are interested in constructing a sample of halo stars, 
we will only consider those stars in the Beers et al. (2000) cat- 
alog with metallicities [Fe/H] < -1.5 dex. 

To minimize any possible contamination from the thick 
disk, we further exclude low metallicity stars near the plane 
and with coplanar and circular orbits, i.e. with disk-like kine- 
matics. Assuming an exponentia l thick disk with scale he ight 
0.8 kp c < h z < 1.5 kpc (Reid & Mai ewskill 19931: iRobin et alJ 
1 19961) and with a velo city ellipsoid (cr g , cr^,, cr z ; ( V rf » = 
(61,58, 39; -36) kms" 1 (TBinnev & Merrifieldlll998l) . we ex- 
clude those stars for which the following conditions are si- 
multaneously satisfied: |Z| < 1.5 kpc, 126 km s _1 < < 
242 kms -1 , \Vr\ < 61 kms -1 , low angular momentum 

L xy = JL\ + lA with respect to L z , namely |L| < 1.1 L z . 



Table 1. Selected sets of metal deficient ([Fe/H] < -1.5 dex) 
halo stars. See text for explanation. 



Type 




1 kpc 


2 kpc 


D 


main-sequence dwarf star 


11 


13 


A 


main-sequence A-type star 


2 


3 


TO 


main-sequence turnoff star 


73 


169 


SG 


subgiant star 


1 


9 


G 


giant star 


78 


133 


AGB 


asymptotic giant branch star 


5 


5 


FHB 


field horizontal-branch star 


5 


5 


RRV 


RR Lyrae variables 


16 


70 


V 


variable star 


1 


3 


TOTAL 




192 


410 



The remaining stars define our selected sample. This sam- 
ple contains 192 stars within 1 kpc and 410 within 2 kpc of 
the Sun. These sets include subdwarfs, giants, and variables, as 
described in Tabled with the following properties: 



- proper motions are available from the HIPPARCOS 
Catalog for 197 stars; the other ground-based measure- 
ments come from the SPM Catalog 2.0, Lick NPM1 
Catalog, STARNET Catalog, ACT Reference Catalog; for 
stars independently measured in two or more catalogs, we 
have averaged the proper motions weighted by their errors; 
in all cases accuracies of a few mas yr~' are achieved. 

- Radial velocities have typical accuracies of the order of 
10 km s _1 . Metal abundances have been determined either 
spectroscopically or from suitable photometric calibrations. 

- Calibrations of absolute magnitude My allow photometric 
parallaxes for which Beers et al. (2000) claim an accuracy 
of ~ 20%. Because of the large distance (d > 200 pc), 
these photometric parallaxes are, formally, more precise 
than the corresponding HIPPARCOS trigonometric paral- 
laxes, however they still remain the main source of uncer- 
tainty in the derived tangential velocities. 

3. Kinematic analysis 

Analytic arguments and high-resolution cosmological simula- 
tions of the formation of dark matter halos (Helmi & White 
1999; Helmi et al. 2002; Moore et al. 2001) suggest that the 
halo should be spatially smooth in the Solar neighborhood (al- 
though see some recent discussion in Diemand et al. 2005; and 
Zhao et al. 2005). This is also supported by measurements of 
the degree of lumpiness in the angular distribution of halo stars 
(e.g. Lemon et al. 2004). However, the same simulations pre- 
dict that while the kinematics of halo objects near the Sun can 
be represented by a smooth multivariate Gaussian distribution, 
the motions of the most energetic particles should be strongly 
clumped and anisotropic. This result has motivated us to ana- 
lyze in more detail the kinematics of the halo samples we iden- 
tified above. 

3.1. Velocity distribution 

The motion of each star in the sample can be specified by the 
velocity vector v = (U, V, W) in 3D linear space, as well as 
in the 2D angular space, by spherical angular coordinates a = 
(cf>,9), where sinfl = VJ\\\ and tan0 = V y /V x . 

Figure shows the kinematic distribution of the selected 
sample within 1 kpc. In Fig. [Ha-c) we plot the U, V, W ve- 
locities. Fig. [fld) shows the distribution of velocity directions 
a = (cp, Q). 

The velocity distribution is relatively smooth, and appears 
to be consistent with a Gaussian Schwarzschild distribution. 
The mean velocities are «[/>, (V +220), <W» = (-26+11,21 ± 
8,-1 + 7) kms 1 and the velocity ellipsoid is radially elon- 
gated, namely (era, cr v , o- w ) = (150 + 8, 106±5,96 + 5)km s" 1 . 
However, a smooth description does not seem to reproduce the 
kinematics of the the fastest objects, as shown by the high- 
lighted points. The 5% fastest moving stars (light color dots) 
seem to be more clumped. The characteristics of the observed 
velocity distribution appear to be in rough agreement with the 
results of CDM simulations as we shall show in Sect. 5. 

We shall perform two different statistical tests on the data 
with the aim of quantifying the presence of large and small 




Fig. 1. Distribution of nearby halo stars in velocity space for our selected sample with [Fe/H] < -1.5 dex within 1 kpc of the 
Sun. (a-c) Velocity projections in the (U, V, W) space; (d) Velocity direction: the position of each particle is given by the spherical 
angular coordinates a = (0, 9) of the velocity vector. Of the 192 stars present in this volume, the 5% fastest are highlighted as 
light color dots. Among them, the asterisks identify the objects/groups with velocity difference less than 42 km s . 



scale anisotropics in the motions of our stars. To establish the 
significance of our results we will compare our results to suit- 
able Monte Carlo simulations. 

Our synthetic data sets have the same number of stars and 
the same spatial distribution as the observed sample. The char- 
acteristic parameters of the multivariate Gaussian used to de- 
scribe the kinematics are obtained by fitting to the observed 
mean values and variances after appropriate convolution with 
observational errors. We generate 100 'observed' samples as 
follows. A velocity is drawn from the underlying multivariate 
Gaussian; it is transformed to a proper motion and radial veloc- 
ity (assuming the observed parallax and position on the sky); 
observational 'errors' of the magnitude described in Sect.2 are 
added to the parallax, the radial velocity and the proper mo- 
tion; these 'observed' quantities are then transformed back to 
an 'observed' velocity. 

The first test we perform on the data consists in quantifying 
the presence of large-scale anisotropics. We implement this test 
by partitioning the 2D angular space ((f), 9) -see Fig.[3d)- into 
cells with roughly similar area. We then count how many stars 
fall in each cell and compare it to the expected number in our 
Monte Carlo simulations. 



Figure |5] shows the results for the partition with 24 cells 
(with an area of 1800- 1500deg 2 ) along4 strips: 0° < \6\ < 45° 
(with = 45°) and 45° < |0| < 90° (with A<p = 90°). 
Plotted are the counts in each cell for the selected halo sam- 
ple within 1 kpc of the Sun (solid line) and for the average 
of 100 Monte Carlo realizations (dashed line). The east-west 
large-scale anisotropy seen in this figure is due to the radi- 
ally elongated velocity ellipsoid in combination with a small 
amount of prograde rotation. 

We find, at the 1 <x level, an excess of stars due to residual 
contamination from the thick-disk (0 ~ 90°) and a lack of stars 
moving toward the SGP (9 90°). To this regard, it is inter- 
esting to note that this deficiency is at odds with some models 
of the Sagittarius dwarf evolution that predict the presence of 
a stellar stream that should cross the Solar neighborhood. The 
stream would be visible as an excess of stars moving toward 
the SGP (e.g. Helmi 2004a; Law et al. 2005). Even though it 
could be argued that this sample may be too bright and possi- 
bly too metal-poor, this result could rule out an oblate or spher- 
ical shape for the Galactic dark-matter halo, and favor a prolate 
mass distribution as suggested by Helmi (2004b). 

In summary, the observed and simulated counts are statisti- 
cally indistinguishable for all partition choices. There is no ev- 
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Fig. 2. ((/>, #)-space number counts for a partition with 24 cells 
with A6 = 45° (mean value 6 C ) and A<p = 45°/90° for the stripes 
along the equator and close the poles respectively. The solid 
line corresponds to the selected sample within 1 kpc of the Sun, 
while the dashed histogram to the average number of counts in 
our Monte Carlo simulations. Error bars are based on Poisson 
statistics. 



idence (at > 2<x-level) for large scale flows crossing the Solar 
neighborhood. 



3.2. Correlation function 

We quantify the deviations from a smooth Gaussian distribu- 
tion due to kinematic substructures by means of the two-point 
correlation function, defined as 



(DP) 
(RR) 



1 



(1) 



where (DD) is the number of pairs of stars in our data with 
velocity difference less than a given value, namely 



{DD) — \ pairs of stars i, j with |v, 



V;|<A 



(2) 



and (RR) is defined analogously for the same number of ran- 
dom points drawn from a multivariate Gaussian distribution de- 
rived from the data set, convolved with expected observational 
errors, and averaged over a hundred realizations. Note that the 
points are not statistically independent, and that our definition 
of £ actually corresponds to the cumulative correlation function 
often used in cosmology in studies of the large-scale structure 
of the Universe (e.g., ICoil et all20ollMullis et all2004l) . 

Based on Poisson counts, we estimate the error of the two- 
point correlation function as 



A* = 



1+j 



(3) 



Thus, in the 3D linear space ( U, V, W) - say 



£(v) measures the excess of pairs of stars moving with veloc- 
ity difference equal or smaller than a given value, above that 
expected from a random sample. Clumping due to kinematic 



Fig. 3. Two-point velocity correlation function for our selected 
sample within 1 kpc (asterisks) and for the 5% fastest moving 
stars (diamonds). In both cases, there is a signal at small veloc- 
ity differences, indicative of the presence of streams. The error 
bars are based on Poisson counts. 



substructures (i.e. groups of stars moving with similar veloci- 
ties) is manifested by an excess at small velocity separations. 

In Figure|3]we show the two-point correlation function £(v) 
for our selected sample within 1 kpc of the Sun (asterisks) and a 
subset comprising the 5% fastest-moving stars (diamonds). We 
use bins of width 21 km s , up to separations of 840 km s -1 . 

The figure shows that, at a 1 cr level, there is a small but 
statistically significant excess of stars with similar velocities 
with respect to what would be expected for a smooth Gaussian 
distribution. The signal in the first bin is weak but significant. 
No correlation is observed at larger separations. 

However the signal is clearly much stronger for the subset 
of 5% fastest moving stars: in this case, the excess of pairs 
of stars with similar velocities is very noticeable, and it is a 
direct indication of the presence of clumps/streams. This signal 
is indeed due to a moving group, which is formed by three stars 
(described in Table|2j and indicated by the asterisks in Fig.^ 

Note also the presence of a certain degree of anti- 
correlation for the interval 200 - 400 km s _1 . This may be due 
to a type of 'clear out effect' as the result of the clumping of 
few objects within the first bins. This effect would have a con- 
siderably smaller amplitude for a larger sample (as is the case 
when our selected sample is considered in its entirety). 

So far we have focused on the 5% fastest moving parti- 
cles within 1 kpc of the Sun. Analysis performed with the 10% 
fastest subsample still shows a significant deviation from a mul- 
tivariate Gaussian, albeit of smaller amplitude. 

Turning to an eight times larger volume (we consider those 
stars within 2 kpc of the Sun) allows us to increase the num- 
ber of stars, but only by a factor of ~ 2. This increase in the 
number of stars, does not necessarily translate into an enlarge- 
ment of the number of streams but could also lead to a better 
representation of each stream. 

However this effect is not obvious in our sample: for the 
larger volume the kinematic group previously singled out loses 
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5 



Name 


[Fe/H] 


a (J2000) 


<5 (J2000) 


D 


U 




V + 220 


W 






(dex) 


(h m s) 


(d m s) 


(kpc) 


(km s~ 


') 


(kms -1 ) 


(km s 


-') 


BPS CS 30339 - 0037 


-2.13 


00 20 28.90 


-36 12 00.7 


0.80 ±0.16 


-292 ± 


54 


73 ±36 


214 ± 


12 


HD 214161 


-2.16 


22 37 08.04 


-40 30 38.4 


0.62 ±0.12 


-290 ± 


23 


123 ± 22 


262 ± 


15 


V* RZ Cep 


-1.77 


22 39 13.05 


+64 51 28.9 


0.41 ± 0.08 


-301 ± 


62 


101 ±27 


240 ± 


47 



one star, BPS CS 30339-0037, which is no longer selected in 
the 5% high velocity tail (simply a consequence of small num- 
ber statistics). Up to velocity difference smaller than 42 km s _1 , 
a weak signal (at the 1 <x level) is noticeable: it is due to the 
two remaining stars of our moving group and a second struc- 
ture with two other stars. Given the low amplitude of this signal 
in comparison to the 1 kpc sample, we are led to believe that 
this second structure is probably a statistical fluke rather than a 
true physical system. 

In summary, we have detected a moving group (with three 
members) among the 5% most energetic stars within 1 kpc of 
the Sun; no other members could be found by increasing the 
volume or relaxing the velocity threshold (from 5 to 10%). 

4. Comparison to the Theory: Simulations 

We now wish to compare our results from the previous sec- 
tion to theoretical models. To this end, we use a high-resolution 
simulation of the formation of a dark matter halo in a ACDM 
cosmology (Springel et al. 2001). From this simulation, we de- 
rive the kinematics of dark matter particles inside spheres of 
2 kpc radius, located at 8 kpc from the Galactic center. In do- 
ing so, we are assuming that these volumes are representative 
of the solar neighborhood. 

In Fig. 0] we plot the kinematics of 2348 particles se- 
lected within one of the 2 kpc spheres. Their velocity distri- 
bution is relatively smooth and appears to be quite consistent 
with a multivariate Gaussian (see Fig. |5}, with principal axes 
(o"i, 0"2, cr 3 ) = (143, 118, 1 10) km s _1 . However, if we focus on 
the motions of the most energetic particles (indicated with gray 
symbols in Fig. [4} this is no longer the case. The 5% fastest 
moving particles are strongly clumped, and their distribution is 
highly anisotropic. 

To quantify the substructures present in this volume we 
compute the velocity correlation function £ (as described in 
Sect. 3). The random comparison sample in this case is the re- 
sult of averaging 100 realizations of a trivariate Gaussian with 
similar moments as the dark matter velocity distribution. In 
Fig. |3 we observe a weak signal in the first bins produced by 
a small excess of particles with similar velocities (asterisks). 
However, if we focus on the 5% fastest moving particles (di- 
amonds), the excess has a much larger amplitude particularly 
at small velocity differences (i.e. Av < 42 km s _1 ), and sim- 
ply reflects the presence of kinematic groups clearly visible in 

Kg. a 

Although the number and the properties of stellar streams 
are likely to be rather different from pure dark-matter streams, 
it is worth noting that these results are qualitatively similar to 



those found for our stellar samples. To try to quantify the de- 
gree of similarity, it is simplest to assume that 10% of the par- 
ticles in these volumes represent stars (i.e. reflecting a "uni- 
versal" baryon fraction). In this way we can randomly define 
"stellar samples" which we subject to statistical analysis, such 
as the velocity correlation function. Such analysis show a weak 
signal at the 1 cr level, which has a smaller amplitude than found 
for the selected stellar sample discussed in Sect. 3. This result 
could suggest that stars are not just a random subset of dark- 
matter particles. This would not be very surprising since stars 
are expected to be much more clustered in the centers of dark- 
matter halos. 

In principle it should be possible to identify which particles 
might represent stars using better motivated physical arguments 
(e.g., by selecting those that have the largest binding energies 
at redshift ~ 10, Moore 2001). However, this is not straightfor- 
ward, as it involves, for example determining the efficiency of 
star formation in each progenitor halo (Robertson et al. 2005). 

In any case, our main limitation in quantifying the degree of 
similarity between the stellar sample and the dark-matter sim- 
ulation lies in the number of particles available in this simula- 
tion. We expect that typ ically, only 1 in 250 particles will be a 
star ( Helmi et al ■1 l 20031) If this is the case, to reproduce our 
selected sample with ~ 400 stars our simulations should have 
10 5 particles in each volume. 

5. Adiabatic invariants 

We are also interested in understanding the properties of our se- 
lected sample in the space of adiabatic invariants, since this is 
generally the preferred space to look for evidence of substruc- 
tures related to past mergers. Here clumping should be strong, 
since all stars originating from the same progenitor have very 
similar integrals of motion, resulting in the superposition of the 
corresponding stream s. In particular, iHelmi et aD 11999) and 
IChiba & Beersl J2000t) have examined the kinematics of metal- 
poor stars in the solar neighborhood and identified a statisti- 
cally significant clumping of stars in the angular momentum 

diagram L z versus L xy = ^L 2 X + Lj. 

The substructure identified by IHelmi et"afl ([l999) con- 
sists of 7 (12) stars with [Fe/H] < -1.6 (-1.0) dex and D < 
1 (2.5) kpc, and was discovered in a sample that was primar- 
ily the result of the combination of HIPPARCOS data with the 
Beers & Sommer-Larsen catalog (1995). Chiba & Beers, us- 
ing the Beers et al. (2000) catalog (i.e. a revised edition of the 

1 This is essentially determined by the ratio of mass in dark-matter 
to that in halo stars for the Solar neighborhood. 
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Fig. 4. Distribution of 2348 dark matter particles from a CDM simulation located in a sphere of 2 kpc radius centered at 8 kpc 
from the Galaxy center, (a-c) Velocity projections in the (vi, Vi, V3) space; (d) Velocity directions. The solid gray circles denote 
the 5% most energetic particles, while the asterisks identify those subsets whose velocity difference is less than 42 km s _1 . The 
kinematics characteristics observed in this set of plots are representative of what is seen in other similar volumes. 
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Fig. 5. Velocity correlation function £(v) for the volume shown 
in Fig. |4] Asterisks and diamonds correspond to £ for the full 
sample and for the subset of 5% fastest moving particles, re- 
spectively. 



previous compilation) confirmed 9 of their stars in the clump 
region which includes also one of the stars in our kinematic 
group, HD214161. However, they later discarded this star for 
having quite different orbital properties than other clump mem- 
bers. 

Here we focus on the space of adiabatic invariants E, L z 
and L xy , although the latter is not fully conserved in an axisym- 
metric potential, and its use may not be generally appropriate 
for the study of substructure in the halo. Nevertheless, we hope 
to obtain further insight into our previous analysis based only 
on kinematics. We assume that the Galactic potential is rep- 
resented by three components Jjohnston et al.lll996l) : a dark 
halo with logarithmic potential, a Miyamoto-Nagai disk, and 
a spherical Hernquist bulge: 



halo 



v;Mn. 



GM dhk 



(a+ 1 



2 +/> 2 ) 2 



(4) 



bulge — 



where M disk = l.O X lO n M , M bu ige = 3.4 x lO lo M ; v h = 
131.5 km s~' and rwo = 300 kpc; a = 6.5 kpc, b = 0.26 kpc, 
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Fig. 6. Distribution of the selected sample of 410 stars with D < 2 kpc in the sp ace of adiabatic in variants. The asterisks denote 
our kinematic group (RHLS), and the grey diamonds the clump identified by Helmi et al. ( 19991 (HWdZZ). These stars are 
used to define a region (limited by the dashed lines in the left panel), which encompasses also other stars with similar angular 
momenta, and which are indicated as diamonds in both panels. CD - 80 328, the candidate from the Nordstrom et al. (2004) 
catalog, has been added and is shown as a triangle. 



c = 0.7 kpc, and d — 12 kpc. This choice of parameters gives a 
circular velocity at the solar radius of 210 km s _1 . 

Figure [6] shows the distribution of the selected sample 
within 2 kpc of the Sun, in the space of adiabatic invariants. 
Note that the sample includes 8 stars (hereafter H WdZZ stars) 
from the stream identified by Helmi et al. ( 1999). These stars 
are shown by gray solid diamonds while our kinematic clump 
(hereafter RHLS) is plotted as asterisks. It is clear that these 
groups have similar momenta (see left panel). If L™ n (L™ n ) and 
L™ ax (L™ x ) denote the minimum and maximum (L z ) values 
for the HWdZZ and RHLS stars, the region Z™ 1 < L xy < I™* 
and Lf m < L- < L™ ax (dotted box) encompasses stars with sim- 
ilar momenta. There are 16 stars located in this region, which 
include five new objects shown as black diamonds in Fig. [6] 
Some of these st ars are also members o f the 'clump-trail' struc- 
ture identified bv lChiba & Beers! feoOO) (hereafter CB). The or- 
bital properties of these 16 stars, based on the Galactic potential 
defined above, are listed in Table [3] The last column serves to 
indicate the membership to the various groups. 

We have looked for additional members of our kine- 
matic group in the Nordstrom et al. (2004) catalog of 
nearby stars. Selecting stars with [Fe/H] < -1.5 dex and 
similar kinematic characteristics as those we identified, we 
find only one possible candidate. CD - 80 328 is a metal 
poor star ([Fe/H] = -1.98 dex) with (U, V + 220, W) = 
(-193, 125, 303) km s _1 and a highly eccentric orbit e = 0.83. 
It is shown as a triangle in Fig. [6] CD - 80 328, having L z - 
1009 kpc km s _1 and L xv = 2422 kpc km s _1 , is located well 
within the box-region defined by the HWdZZ stars. However, 
given its much lower binding energy (E = -0.84- 10 5 km 2 s~ 2 ), 
it is more likely that CD - 80 328 is a member of our moving 
group. 

According to numerical simulations carried out by 
iBrook et all i 20031) . stars from an accreted satellite show typ- 
ically highly eccentric orbits (e > 0.8) and strongly correlated 
velocities. In view of this, the stars of our kinematic group 



(with 0.81 < e < 0.85) could well share a common origin, 
and be stellar debris from an accreted satellite. 

It is worth noting that, while having similar momenta, our 
kinematic structure and the HWdZZ clump have somewhat dif- 
ferent energies. This could imply that these groups have differ- 
ent origin. However, it seems also plausible that these groups 
have been stripped off from the same progenitor at different 
times. In this scenario, HWdZZ stars should have been released 
in a later galactic passage than those in our kinematic clump. 
The difference in orbital energy between the groups could be 
the result of different binding energies of the progenitor galaxy 
(higher and lower, respectively) due to the effects of dynamical 
friction upon this system while orbiting the Milky Way. 

6. Conclusions 

We have extracted samples of metal -poor ( [Fe/Hl < -1 . 5 dex) 
halo stars in the Solar neighborhood from Bee rs et alJ (|2000) 
catalog, and carried out a statistical analysis of their kinemat- 
ics. 

Based on clustering in the velocity space, we have found 
evidence of substructures in the motions of the fastest moving 
stars, at a level which seems to be consistent with that predicted 
by high resolution simulations of dark matter halos in a hierar- 
chical universe. 

The moving group responsible for this signal is comprised 
by three stars, whose kinematic and metallicity characteristics 
are similar to the streams found by Helmi et al. (1999), albeit 
on somewhat more loosely bound orbits. 

Our sample of halo and high velocity stars is too small to 
make definite statements about the importance of accretion in 
the formation of the Galactic halo. This would require a sam- 
ple with a few thousands nearby halo stars (i.e., ~ 10 times 
larger than our selected sample) with accurate space velocities. 
Such sample sizes will become available in the near future, 
thanks to spectroscopic surveys like RAVE (Steinmetz 2003) 
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Table 3. Characteristics of the stars located in the region L™ n < L„. < L™ ax and L™ n < L z < L™ x defined by the Helmi et al. 
(1999) structure (within the box shown in the left panel of Fig.|6]) 



Name 


[Fe/H] 


E/10 5 


U 


Lxy 


e 


Membership 




(dex) 


(km 2 s- 2 ) 


(kpc km s -1 ) 


(kpc km s -1 ) 






BPS CS 22948-0093 


-3.72 


-0.57 


687 


1805 


0.92 




HD 214161 


-2.16 


-0.73 


952 


2173 


0.81 


RHLS 


V* RZ Cep 


-1.77 


-0.77 


710 


1943 


0.85 


RHLS 


BPS CS 30339-0037 


-2.13 


-0.87 


586 


1909 


0.81 


RHLS 


BPS CS 22189-0007 


-2.12 


-0.97 


683 


1624 


0.77 




BPS CS 29513-0031 


-2.79 


-1.00 


1171 


2189 


0.47 


CB 


BD+10 2495 


-1.83 


-1.02 


1318 


2140 


0.47 


HWdZZ, CB 


HD 119516 


-2.49 


-1.04 


1212 


2055 


0.50 


HWdZZ, CB 


V* TT Cnc 


-1.57 


-1.06 


562 


2483 


0.42 


HWdZZ 


CD-36 1052 


-2.19 


-1.06 


1078 


2132 


0.45 


HWdZZ, CB 


HD 237846 


-2.63 


-1.09 


774 


1443 


0.71 


HWdZZ 


BPS CS 29504-0044 


-2.04 


-1.10 


958 


2091 


0.42 


CB 


V* AR Ser 


-1.78 


-1.12 


322 


1934 


0.46 


HWdZZ 


V* TT Lyn 


-1.56 


-1.11 


748 


2124 


0.45 


HWdZZ, CB 


HD 128279 


-2.20 


-1.18 


1194 


1844 


0.23 


HWdZZ, CB 


BPS CS 22876-0040 


-2.20 


-1.19 


922 


1561 


0.42 


CB 



Membership: RHLS: Re Fiorentin et al. (this paper); HWdZZ: lHelmi et alJll999li : CB: IChit 



and SDSS-II/SEGUE (Beers et al. 2004), which could be com- 
bined with proper motion catalogs such as UCAC2 (Zacharias 
et al. 2004), SPM (Platais et al. 1998), GSC-II (McLean et al. 
2000), USNO-B (Monet et al. 2003), to obtain full phase-space 
information. 

The existence of structures in the halo, if confirmed by fur- 
ther studies, is of great importance for constraining models of 
the formation and evolution of the Galaxy. The space astromet- 
ric mission Gaia (Perryman et al. 2001) will collect samples 
of millions of stars in our Galaxy as well as in our nearest 
neighbours with very accurate positions, proper motions, and 
trigonometric parallaxes which will dramatically improve the 
situation, and revolutionize our knowledge of the Galaxy. 
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